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Abstract 


Perturbative QCD predicts the behavior of scattering at high energies and 
fixed (sufficiently large) transferred momenta in terms of the BFKL pomeron 
(or short distance pomeron). We study the prospects for testing these predic¬ 
tions in two-photon processes at LEP200 and a possible future e ± e _ collider. 
We argue that the total cross section for scattering two photons sufficiently 
far off shell provides a clean probe of BFKL dynamics. The photons act as 
color dipoles with small transverse size, so that the QCD interactions can be 
treated perturbatively. We analyze the properties of the QCD result and the 
possibility of testing them experimentally. We give an estimate of the rates 
expected and discuss the uncertainties of these results associated with the 
accuracy of the present theoretical calculations. 
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I. INTRODUCTION 


The behavior of scattering in the limit of high energy and fixed momentum transfer 
is described in QCD, at least for situations in which perturbation theory applies, by the 
Balitskii-Fadin-Kuraev-Lipatov (BFKL) pomeron |l[ (or short distance pomeron). Attempts 
to test experimentally this sector of QCD have started in the last few years, mainly based 
on measurements of deeply inelastic events at low values of the Bjorken variable x in lepton- 
hadron scattering J2|] and jet production at large rapidity separations in hadron-hadron 
collisions In this paper we study the possibilities for investigating QCD pomeron effects 
in a different context, namely in photon-photon scattering at e + e~ colliders, where the 
photons are produced from the lepton beams by bremsstrahlung. The results of this paper 
also apply to e~ e~ or /i ± /U colliders. Some aspects of this study have been presented in 
Refs. Q| [n . 

The quantity we consider is the total cross section for off-shell photon scattering at 
high energy. This can be measured in e + e~ collisions in which both the outgoing leptons 
are tagged. This cross section presents some theoretical advantages as a probe of QCD 
pomeron dynamics compared to the structure functions for deeply inelastic scattering off a 
proton (see for instance the discussion in Ref. 0) or a (quasi)-real photon (see for example 
Ref. [7]), essentially because it does not involve a non-perturbative target. Unlike protons 
or quasi-real photons, virtual photon states can be described through perturbative wave 
functions. In some respects the off-shell photon cross section presents analogies with the 
process of scattering of two quarkonia (or “onia”), which has been proposed as a gedanken 
experiment to investigate the high energy regime in QCD ||. In this case, non-perturbative 
effects are suppressed by the smallness of the onium radius. In the case of virtual photons 
the size of the wave function is controlled by the photon virtuality instead of the heavy 
quark mass. It is an interesting feature of investigations at e + e _ colliders that this size can 
be tuned by measuring the momenta of the outgoing leptons. 

On the other hand, such experimental studies may prove to be difficult due to the 
smallness of the available rates. As we shall see, for large photon virtualities Q 2 the cross 
section falls off like 1 /Q 2 . An estimate of the number of events that one may expect to be 
available for such studies at LEP200 and a future linear e + e~ collider will be provided later 
in the paper. 

There have recently been other investigations of the high energy regime in the context 
of photon-photon scattering. Balitskii [j|] has proposed an expansion of the scattering am¬ 
plitude in the high energy limit in terms of Wilson line operators. This method provides 
an elegant reformulation of the BFKL problem, and may prove to be useful to get beyond 
the leading logarithm approximation. Bartels, De Roeck and Lotter [|Tf| have evaluated the 
photon-photon cross section at present and future e + e _ colliders. Their results are similar 
to those in Refs. [|] []|. In Ref. [7j] detailed studies have been carried out for diffractive meson 
production and photon structure functions at LEP200. 

To describe the electron-positron scattering process (Fig. [I]), we will parametrize the 
five-fold differential cross section as 
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( 1 . 1 ) 


da^ e+e ) 

dxA dx B dQ\ dQ 2 B d(f >/(2 7 r) 


d cr( e+e ) 

dxA dx B dQ\ dQ 2 B 


1 + Ai cos (p + A 2 (2 cos 2 0 — 1) 


Here we denote by xa, x b the fractions of the longitudinal momenta of the leptons A and 
B that are carried by the photons, by Q\ = Q\ — ~Qb ( 1bh the photon virtualities, 

and by 0 the angle between the lepton scattering planes in a frame in which the photons 
are aligned with the z axis. The overall factor in the right hand side gives the distribution 
averaged over the angle 0, while Ai, A 2 are the asymmetries. 



FIG. 1. The photon-photon scattering process in e + e collisions. 


We will start by considering the 0-averaged cross section, and we will express it in 
the equivalent photon approximation [|ll|] by folding the 7* 7* cross section with the flux 
of photons from each lepton. This flux is proportional to the probability density for the 
splitting e —> ey, and depends on the photon polarization. We will thus write 


Qa Qb da^ e+e } 
dxA dx B dQ\ d,Q 2 B 



{x B ) 4*7* (xa x b s,Q 2 a ,Q b ) 


+4 T /l+ ( x A ) P^l- (: x B ) cr^) (x A x B s, Q% Q\ ) 

+ P 7/'e ) + ( X A) P^/l- (X B ) <T§ 7* (XA X B S, Q\, Q%) 

+ P§1+ {xa) P§1- {x B ) 4*7* { X A x b S, Qa, Qb) } 


( 1 . 2 ) 


where the transverse and longitudinal photon flux factors P (T 0 are given by 



1 + (1 — x) 2 
x 



1 — x 

x 
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(1.3) 


and 44* > with a,b = T,L, is the cross section for the scattering of two photons with 
polarizations a and b. 

We will proceed in the following way. In Sec. || we set our notations and describe the 
Born approximation to the photon-photon cross section at high energy. In this section we 
concentrate on the case of transversely polarized photons, that is, the cross section 44* 
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Eq. ( |1.2|) above. In Sec. jnljwe extend these results to include the full polarization dependence 
and we discuss the associated asymmetries. In Sec. |T^ we consider the summation of the 
leading logarithmic corrections to the photon-photon cross section due to BFKL pomeron 
exchange, and emphasize how the perturbative results depend on the total energy and the 
photon virtualities. Secs. |V| and 0 are devoted to discussing some of the limitations of the 
treatment based on the BFKL equation. In Sec. 0 we focus on the dependence of the cross 
section on two mass scales (in the running coupling and in the high energy logarithms), which 
are left undetermined in a leading log analysis. In Sec. 0 we consider the limitations of 
using the BFKL approach that follow from the behavior at very large s. Sec. 01 illustrates 
how to relate the summed result for the photon-photon cross section to the small-x behavior 
of the photon deeply inelastic structure function. In Sec. |VIII| we compare the QCD result 
with expectations based on traditional Regge theory. In Sec. [IX] we consider the limit of low 
photon virtualities and discuss the region of transition between hard and soft scattering. The 
rates at the level of the e + e _ cross section are examined in Sec. [X]. Some concluding remarks 
are given in Sec. |Xj. We collect in Appendix |X] the details of the Born order calculation, 
and in Appendix |B] some formulas which are useful for comparing the gluon-exchange and 
quark-exchange contributions to the photon-photon cross section. 


II. NOTATIONS AND LOWEST ORDER CALCULATION 

We consider the total cross section for the scattering of two transversely polarized virtual 
(space-like) photons 7 *(< 74 ) and 7 *(ge), with virtualities q\ = — Q\ and q 2 B = — Q B , in the 
high energy region where the center-of-mass energy -y/s = J ( q a + qs ) 2 is much larger than 
Qai Qb • We also suppose that the photon virtualities are in turn large with respect to 
the QCD scale A q CD , so that the process occurs at short distances (much smaller than 
A qcd ~ 1 bn) and QCD perturbation theory applies. 

We work in a reference frame in which the incoming photons have zero transverse mo¬ 
menta, and are boosted along the positive and negative light-cone directions. For a four- 
momentum p M , we define the “+” and ” momentum components p + , p~ as 

p ± = (p°±p 3 )/V2 . ( 2 . 1 ) 

The incoming photon momenta are parametrized as follows in a notation where q M = 

(<? + , q~, qr)= 

q “ A= ■ (2 - 2) 

Here the “+” and ” components q\, q B approximately build up the total energy s, and 
are much larger than the initial virtualities 

~ s • ( 2 . 3 ) 

The Born contribution to this cross section is given by the high energy approximation 
to the reaction 
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7 *(u) + 7 *(Qb) -> q(pA) + ?(pa) + q(Pb) + g(p. 


B 


(2.4) 


in lowest order perturbation theory. The corresponding diagrams involve the exchange of 
two gluons between the two quark-antiquark pairs |12|, and are exemplified in Fig. 

We parametrize the outgoing quark momenta as 


Pa = (j'Aq'Ai z'aQbi Pa 


Pb = {z'r ( A^ z b ( 1biPb 


(2.5) 


and denote by the exchanged gluon momentum. 

In the high energy limit defined by Eq. (2.3), the kinematic region which dominates 
the integrations is the one in which the transverse momenta flowing in the loops are of the 
order of the initial virtualities, and the light-cone components of the exchanged gluon are 
suppressed with respect to the transverse momenta by a quantity of order Q A /y/s or Q B /^/s, 
so that k 2 ~ —k 2 . 



FIG. 2 . One of the two-gluon exchange graphs contributing to the high energy 7 * 7 * cross 
section in the Born approximation. The gluons can be attached to the quark lines in 2 4 different 
ways. 


The squared amplitude for the graph in Fig. [2j. integrated over the final quark and 
antiquark phase space, for fixed transverse photon polarizations £a, £b, is given by 




d 4 k d 4 pA d 4 p B 1 


(2 7r ) 4 (2 7r ) 4 (2 7r ) 4 {k 2 Y 
x 2 7id + ((q A - Pa~ k) 2 ^) 2n5 + ((q B - p B + k) 2 ^j 


2?r <5 +(pa) 27r <5+(pi) 

a,b 


( 2 . 6 ) 


TV 


x 


i>A (* 9st r 7c) i (jf>A + ft) (* e a e j A ) {f A - 


-ig s t s lp){-i)^A ~ 4 A )(- ie aeft A ) 


Tr 


x 


((Pa + k) 2 + i e) ((p A - q A ) 2 - i e) 

ft> B ( i 9st r 7 ") i C i> B - ft) (ie b ef B ) {j B ~i> B + ft)(- i 9st s -f 0 )(-i){]/) B - 4 B )(-ie b <4 B ) 


((; Pb ~k) 2 + i e) ((p B - q B ) 2 - i e) 


Here e a , are the electric charges of the quarks in units of e = V4 7r a, and the indices a, b 
run over the light quark flavors, u, d, s. (The flavors c, b need a separate treatment.) In the 
high energy approximation this amplitude takes the form 
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C dzA l' dzB wr 


\M\ 2 — 32a 2 a 2 fee?) ( 2 ^( 2 *) / 


d 2 k d 2 pA d 2 ps 

(2^F 
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x 


, 7 . (2tt) 2 (2vr) 2 

[4:Z A (1 - z A )s A • p 4 £a • (pa + k) - Pa ■ (pa + k)] 

- z A ) Q a + p\] 

P b ■ (Pb - k)] 


[za(1- z a )Qa 
[4z.b(1 — zb)sb 


A + (P-4 + k) 


[za( 1 
Ps £b ■ (Pb - k) - 


[z B (1 - z B ) Qb + (pb 


k ) 2 ] [zb{1~ zb)Qb + Pb] 


(2.7) 


See Appendix A for details of this calculation. 

The amplitudes for the other graphs, in which the gluons are connected to different 
fermion lines, can be derived from Eq. (|2.7|) by using the replacements 


za(1-za)Qa+ (Pa + k) 2 


f -Z A 

ZA 


Z A (1 - Z A ) Q\ + p 2 4 


( 2 . 8 ) 


in the denominator, and 

(Pa + k) -> - 1 - p A (2.9) 

in the numerator. (Analogous replacements hold for the momentum components of the 
quark p B ). In addition, the contribution from the graphs in which the quark and antiquark 
lines are interchanged can be obtained by symmetrizing the above expressions with respect 
to pa —> Pa + k (and, analogously, p B —»► p B — k). 

We add the graphs, and divide by 2 s to form the 7* 7* cross section 

77. = 2_|X| 2 . (2.10) 

We find 

^.(s,Ql,Q 2 B ,e A ,e B ) = 2-J ^-^-G(k',Q 2 A ,e A )G(-h,Ql,e B ) . (2.11) 

The factors 1/ (k 2 )" in this formula come from the gluon propagators. These factors multi¬ 
ply functions G(k; Q 2 ,e), which describe the coupling of the exchanged gluon to the quark- 
antiquark system created by the virtual photon with virtuality Q 2 and transverse polariza¬ 
tion e. The explicit expression for G is 


G(k; Q 2 ,e) = 4a a s ^e 2 j J ~~ J Q dz (2.12) 

[ [p 2 ~ 4z(l — z)(g • p) 2 ] [p • (p + k) - 4z(l -z)e-pe-(p + k)] ) 

1 [p 2 + z (1 - z) Q 2 } 2 [(p + k) 2 + z{l^*,z) Q 2 } [p 2 + z(l — z) Q 2 } J 


The functions G can be thought of as “color functions” of the virtual photon since they 
describe the color flow in the qq states. From the point of view of light-cone perturbation 
theory 0 , they correspond to the coupling of the null-plane photon wave function to gluons. 

In the remainder of this section we discuss the case of the average over the two transverse 
photon polarizations. The detailed polarization dependence of the color functions G and 
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the associated polarization asymmetry in the cross section are treated in Sec. 0- By taking 
the polarization average 


1 

2 


£ 4 a ) 4 a> 



(2.13) 


we define the function G i 

G 1 (k 2 /Q 2 ) = i^ G ( k i<3 2 . £|A) ) • (2.14) 

2 A 


We find 


G 1 (k 2 /Q 2 ) =4aa s 


d 2 p 


d z 


i r Jo 


z 2 + ( 1 - zf 


X 


p z 


p ■ (p + k) 


[p 2 + z(l- z) Q 2 Y [(p + k) 2 + z(l- z) Q 2 } [p 2 + z(l"-z) Q 2 } 


(2.15) 


Each one of the two p-dependent terms in the integrand of Eq. ( [2.1 5|) would lead to an ultra¬ 
violet divergent integral, but the divergence cancels in the sum, illustrating that the gluon 
does not couple to the color singlet qq system in the limit p 2 —> oo. The photon virtuality 
Q 2 regularizes the denominators in the infrared region, p 2 —> 0. In the limit of small Q 2 , the 
probability density Pq/^z) = [z 2 + (1 — z) 2 ] /2 for the splitting of a transversely polarized 
photon into a quark-antiquark pair (7 —> qq) correctly factors out in front of the logarithmic 
singularity dp 2 /p 2 associated with the region of strong ordering Q 2 ■ C p 2 -C k 2 . 

By introducing the Feynman parametrization 


[(p + k) 2 + z (1 - 


1 

- z) Q 2 ] [p 2 + 2 (1 - z) Q 2 } 1 

_ 2 (1-A)_ 

[(p + A k) 2 + 2 (1 - z) Q 2 + A (1 - A)k 2 ] 2 


(2.16) 


and carrying out the integration over the shifted transverse momentum variable p ; = p +A k 
in Eq. ( ^.15| ), one can obtain the following useful representation of the function Gi(k 2 /Q 2 ) 
as an integral over two dimensionless variables: 


Gd(k 2 /Q 2 ) 



k 2 



[A 2 + (l-A) 2 ] [^ 2 + (l-^) 2 ] 
A (1 — A) k 2 + ^ (1 — 2 :) Q 2 


(2.17) 


This representation explicitly shows that the distribution (k 2 ) 1 Gi(k 2 /Q 2 ) is symmetric 
under interchange of the space-like boson virtualities k 2 and Q 2 . I 11 the configurations 

in which one of the virtualities is much smaller than the other one, this distribution is 
logarithmically enhanced. More precisely, we find 


— Gi(k 2 /Q 2 


4 a a. 



2 1 
3 g 2 


o 2 

ln^ + 0(l) 


k 2 < Q 2 , (2.18) 
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k 2 


Gi(k 2 /Q 2 ) ~4a« s 



h 4 +o{i) 


k 2 »Q 2 


(2.19) 


where the coefficient in front of the logarithm is given by the first moment of the splitting 
density 


d x 


x 2 + (1 — xY 


2 

3 


( 2 . 20 ) 


The logarithmic behavior at small k 2 comes from the region k 2 -C p 2 -C Q 2 . Here the 
quark transverse momentum is much smaller than the photon virtuality, and the quark 
longitudinal momentum fraction is very small, z -C 1 (or (1 — z) <C 1). This region is 
sometimes referred to as the aligned-jet region, and corresponds to configurations in which 
the qq system fluctuates to large sizes [[TjJ. Note that while for the case at hand of the 7* 7* 
total cross section this region contributes only a logarithmic enhancement, for the case of 
non-inclusive processes, such as processes involving rapidity gaps, this is expected to become 
the dominant contribution |TS| . 


The integration over the parameter A in Eq. (|2.17|) can be explicitly performed. This 
yields 


Gi{rj) = 2aa s 



0 a/1 - ^ 


1 - 


C 


V + £/2 ]n / Vh + £ + yV 
_yjv (V + 0 n Wv + Z-y/ri, 


- 1 


where we have introduced the variables 
£ = 4z{l-z) 


V = k 2 /Q 2 


( 2 . 21 ) 


( 2 . 22 ) 


In Fig. we report the result of the numerical evaluation of the integral ( p.21| ) by plotting 
the function G\ versus rj = k 2 /Q 2 . 

The unpolarized 7* 7* cross section to Born order, , can be obtained by inserting 


the result for the function G\ in the general formula (|2.11|) . It is convenient to use the 
representation (|2.17|) of G±. By substituting this representation in Eq. (|2.11|) , and carrying 
out the integration over the gluon transverse momentum k, we obtain 


cr 


*(s, Q\, Q 2 b ) — 


2 a 2 a 2 (Eye 2 )' 


7r 


to 

-1 


rl 


dz A 

Jo 

[4 + (1 - -a) 2 ] 

_ 0 

[a^ + (1 - Aa) 2 ] 


x 


dz B z 2 B + (1- z B f 


dX r 


A| + (1 - A s ) s 


x 


In [Q\ z A (1 ~ z A ) X B (1 - A b)/ {Q\ z b (1 - z B ) A^ (1 - Aa))] 
Qa z a (1 — za) X b (1 — \ B ) — Q 2 b z B (1 — z B ) Aa (1 — Aa) 


(2.23) 


This formula provides an expression for the cross section in the large s limit in terms of 
dimensionless integrals, which allows us to study the dependence on the energy and mass 
scales. To this order in perturbation theory the cross section has a constant behavior with 
the energy s. The cross section depends on the mass scales Q\ and Q 2 B only. Factoring out 
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k 2 /Q 2 

FIG. 3. The k 2 /Q 2 dependence of the color functions G\ (see Sec. ID and G *2 (see Sec. |I|). 
We plot the normalized functions G = G/(4aa s J2 q e l)- 



an overall scale factor 1 /(QaQb) in Eq. (]2.23|) , we are left with a function a (r) of the ratio 
r = Qa/Qb- 


a 


(o) 



°( r ) 

QaQb 


(2.24) 


The function a(r) can be computed by performing numerically the integrations over the 
dimensionless variables z’s and A’s. The result is plotted in Fig. [1|. 

Notice that the dependence of a(r ) on r is rather mild when r is near 1. Thus, in this 
region, a(r) in Eq. ( |2.24| ) could be treated as a constant, so that 

a (0) ~l /(Q a Q b ) , r ~ 1 (2.25) 


When, in contrast, the two virtualities are widely disparate from each other, either r> 1 
or r < 1, the function a(r ) has a strong dependence on r, vanishing linearly (modulo 
logarithmic enhancements) with either r -1 or r, respectively. This provides the overall 
change of scale in the cross section 


cr (0) ~ 1/Qa , 

q 2 a ^>qI , 

(2.26) 

n (0) ~ 1/Qb , 

Qb > Qa ■ 

(2.27) 


We will come back to this later on and examine the details of these behaviors numerically 
(see Sec. |V'11| ). 


9 











QJ Q b 

FIG. 4. The high-energy 7 * 7 * cross section in the Born approximation as a function of the 
ratio r = Qa/Qb between the photon virtualities. The solid line is the rescaled cross section a(r) 
defined in Eq. (2.24). The dashed line is the polarization asymmetry discussed in Sec. H (see 
Eq. ( |3 .1 2|) ), multiplied by 10 2 . 


III. POLARIZATION DEPENDENCE 


We now study the dependence of the color functions G and of the photon-photon cross 
section on the polarization of the virtual photon. Working in the frame defined by Eq. 
and denoting by u , v the lightlike unit vectors 

u = (1,0,0) , v = (0,1,0) , (3.1) 


we introduce the following decomposition of the polarization tensor for the photon qa- 


V- if N if N Qa W + v fl q Ah 

9HV = 2^ £ M A ) £ ^9A) - T- —VnV* 


*= 1,2 


(Qa-v 


2 


q. a • v 


(3.2) 


The index i in the Erst term in the right hand side of Eq. (|3.2|) runs over the transverse 
polarizations. We describe these polarizations by using the linear basis 


= K ,i = 1,2 . 


(3.3) 


From the second term in the right hand side of Eq. (|3.2|) , we define the longitudinal polar¬ 
ization vector as 


Lf \ Qa 
£ Ma) = - V„ 


qA ■ v 


(3.4) 
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The last term in Eq. ( fh2| ) does not contribute because of current conservation. We thus 
have 


~ 9^ = J2 £ ]m{Qa) £ U r u) ~ e%(q A ) £ t(qA) + {gauge terms} . (3.5) 

i=1,2 

Formulas analogous to Eqs. ( |T2l) - (|3.5|) , but with v M replaced by u hold for the photon q B . 

Let us first consider the case of transverse polarizations. It is convenient to introduce 
the tensor amplitude Q ltv by rewriting Eq. (|2.12|) in the form 

G(k ] Q 2 ,e) = g^(k ] Q 2 ) £fl e u . (3.6) 


The explicit expression for Q can be found in Appendix A. For polarization indices a, a' = 
1,2, L we also introduce the notation 


g 


E a E a 


G L 


(3.7) 


We can parametrize the transverse components of the tensor G aa ' in terms of two scalar 
functions G\ and G 2 \ 


G H ' (k; Q 2 ) =5™' Gi(k 2 /Q 2 ) 




G 2 (k 2 /Q 2 ) 


i,i' = 1,2 . (3.8) 


The function G± represents the unpolarized color function, which we have discussed in the 
previous section (Eqs. ( [2.14 ), ([2. 15 )). The function G 2 carries the information on the 
polarization dependence. Its explicit expression is 


G 2 {k 2 /Q 2 ) = 4aa s J dzAz(l — z) 

p 2 - 2 (p ■ k) 2 /k 2 P ' (p + k) — 2 (p + k) • k p ■ k/k 2 1 

[p 2 + z{l-z) Q 2 } 2 [(p + k) 2 + z(l — z) Q 2 ] [p 2 -\- z{l — z) Q 2 ] \ 

The integrals in Eq. (|3.9|) can be handled using the same procedure as in the unpolarized 
case: 


G 2 (k 2 /Q 2 ) = 16 aa s I 

V q 

= 2 aa s 55 C 


fi 


r* i 


A (1 — A) z (1 — z) 


k 2 I d z I d A , . ,,, „ . 

Jo Jo A (1 — A) k 2 + 2 (1 — z) Q 2 


(3.10) 


"1 


dt 




e/2 


' Vv + t + Vv " 
yJr] (77 + 0 \Vv + Z- Vv, 


In 


where we have used the variables defined in Eq. (|2.22|) . 

Correspondingly, the cross section for scattering transversely polarized photons, 
Eq.( [2.11|) , can be decomposed as 


a, 


77 


(s, Q 2 A, Qb, £ A, £b) = cr(s, Qa , Q 2 b) [l + (2 {e A ■ £ b) 2 ~ l) A(s, Q% Q 2 B )\ . (3.11) 


The polarization average a has been given to order a 2 in the previous section (Eq. (|2.23|) ), 
and the asymmetry A to the same order a 2 is given in terms of the color function G 2 as 
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«4( s ) Q\, Q 2 b ) — 


d 2 k 1 


a(s,Q 2 A ,Q 2 B ) 16 tt J 7T (k 2 ) : 


■G 2 (k 2 /Q^)G 2 (k 2 /Q|) 


(3.12) 


Numerical results are reported in Figs. |3| and [|, where we plot G 2 versus r/ = k 2 /Q 2 , 
and show the dependence of A on the incoming photon virtualities. Unlike G\ 1 G 2 is not 
logarithmically enhanced in the regions k 2 < Q 2 , k 2 Q 2 . This is related to the fact that 
the splitting function associated with G 2 has zeros at the endpoints of the spectrum in the 
longitudinal momentum fraction, z = 0 and z = 1 (see Eq. ( 3.1 0|) ), whereas the unpolarized 
splitting function goes to a finite constant. The asymmetry A associated with the photon- 
photon scattering process, calculated here in the Born approximation, contributes to the 
asymmetry A 2 in Eq. ( |1.1| ) at the level of the e + scattering. Numerical results for the 
e + e~ process will be given in Sec. [X]. 

We now move on to the case of the longitudinal polarization. Let us consider the 
longitudinal-longitudinal color function: 


G ll = g 


H v 




(3.13) 


Following the lines of the calculation described in detail for the case of transverse photons, 
we find 


G 


x 


LL 


(k 2 /Q 2 ) = 4cta s ^e 2 j / ~~ j 0 dzAz 2 (l-zf 


Q 2 


Q 2 


[p 2 + z (1 - z) Q 2 Y [(p + k ) 2 + 2 (1 - z) Q 2 } [p 2 + z(l-z) Q 2 } 


(3.14) 


This contribution equals the color function G 2 given above. This can be seen by introducing 
an integral over a Feynman parameter A, then integrating over the transverse momentum p 
in Eq. (|3.14| ). We get 


G LL (k 2 /Q 2 ) = 16 a a s e 2 \ k 2 f dz f d X 


A (1 — A) z (1 — z) 
o A (1 — A) k 2 + ^ (1 — ^) Q 2 


= G 2 


(3.15) 


As noted above, G 2 has no logarithms at small k 2 . This corresponds to the absence of 
aligned-jet terms for longitudinally polarized photons ]15|. Contributions from longitudinally 
polarized photons enter the e + e~ cross section ( 0 ). They will be included in the numerical 
estimates that we give in Sec. [X|. 

Finally, we consider the interference contribution between longitudinal and transverse 
polarizations: 


G iL =G Li = g^ £ ^ £ l ; i = 1,2 . (3.16) 

In the high energy approximation in which we are working, this contribution vanishes. This 
can be seen explicitly by writing the color function in the general form 

G iL = G Li = G 3 (k 2 /Q 2 ) k'7|k| , (3.17) 

and computing the invariant function G 3 . We get 
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G 3 (k 2 /Q 2 ) = 16 aa s 



|k| Q dz d A 

Jo Jo 


\(l-\)z(l-z)(l-2z) 

A(l-A)k 2 + z(l-^) g 2 


(3.18) 


Interference terms of the kind in Eq. ( p. 17] ) would give rise to the asymmetry A\ (see 
Eq. (|1.1|) ) at the level of the e + e~ scattering process. We thus see that this asymmetry 
vanishes in the high energy approximation. 


IV. SUMMATION OF LEADING LOGARITHMS 


We have seen that the two-gluon exchange mechanism gives rise to a constant 7 * 7 * total 
cross section at large s, cT 0 )(s, Q\, Q 2 B ) ~ a 2 a 2 f(Q \, Q B ). To higher orders in perturbation 
theory, the iteration of gluon ladders in the f-channel promotes this constant to logarithms, 
and the perturbative expansion of the cross section at high energy has the form 


a. 




~ a 


(o) 


OO 

1 + flfc (<U L ) + 

k =1 


L = ln(s/ Q 2 ) 


(4.1) 


where Q 2 is a scale of the order of the initial photon virtualities, the sum represents the 
series of the leading logarithms to all orders in the strong coupling a s , and the dots stand 
for non-leading terms. 

To study the high energy behavior, it is convenient to analyze the cross section in its 
Mellin-Fourier moments, defined by 


/ D 2 n 2 \ — [ a+l °° d^ NL ir\2 n 2 \ 

J a—ioo Z71X 


(4.2) 


where the iV-integral goes along a contour parallel to the imaginary axis and to the right 
of any singularities in ctjv- We see from this definition that a constant behavior of the cross 
section with the energy s is generated by a simple pole in the moments ctn at N = 0, while 


powers of logarithms are generated by multiple poles at N = 0. The inverse of (4.2) is 


°n(q 2 a ,qI) = * {if) 


(4.3) 


To sum the leading logarithmic terms, the basic observation of BFKL [l| is that the 
logarithms arise when multiple soft gluons are emitted into the final state. These gluons 
have transverse momenta k± of the same order as Qa and Qb and have strongly ordered 
rapidities, lying between the rapidities of the quarks in photon A and the quarks in photon 
B. Along with emission of real gluons, one also includes the exchange of corresponding 
virtual gluons. 

Fig. j| illustrates how the two gluon exchange graph of Fig. [2| is generalized to allow for 
multiple gluon emission. The quarks comprising photon A couple to a gluon with momentum 
k 1 ^, while the quarks comprising photon B couple to a gluon with momentum k B . Gluons 
can be exchanged or emitted into the final state inside the subgraph labeled T. In fact, 
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the gluons that carry momenta k^ and kg are, in general, combinations of soft gluons that 
carry a net color octet charge, that is, reggeized gluons [l|. From a kinematic viewpoint, 
we can treat these as being equivalent to ordinary perturbative gluons. We consider the 
unpolarized 7 * 7 * cross section and adopt the following notation for the diagram in Fig. ||: 



dk\ J{-k\/q\,-^A\Q 2 A ) 


X 


r d 2 k B riB 

J 7rk^ Jo 


dk B J(k B /q B ,]tB' 1 Q 2 B ) G(k~A, k B ,\s.A,k-B) ■ 


(4.4) 


Here Q represents the four-point Green function for gluons k A and k B , while the functions 
J represent the quark loops. We have made the following approximation. We note that the 
quark loop for photon A depends sensitively on the minus component, of the momentum 
that enters the quark loop via gluon A. On the other hand, the quarks in photon B have 
very large minus components of momenta. Furthermore, because of strong rapidity ordering, 
all the other gluons inside of T have minus components of momenta that are much larger 
than k Thus we can neglect kJ everywhere except in the quark loop for photon A. Then 
we include the integration over k^ in the definition of J(— k\/q\, — k 4 ; Qa)- Similarly we 
neglect k~J, everywhere except in the quark loop for photon B and we include the integration 
over k B in the definition of J[k B / q B , kg; Q B ). 



FIG. 5. Factorized structure of the virtual photon cross section in the high energy limit. 

Consider Eq. Qf.4|) specialized to the case of two gluon exchange. The approximate gluon 
four-point function Q is trivial in this case: 

Qo(k\i k B , k A , k B ) = ^ <5(k 4 - k B ) 5{k\) 5{k B ) . (4.5) 

We identify 


J (0, — k^; Q\) = G, (k^/ Q\) 
Then we recover the earlier result (|2.11|). 


J(0,k B ;Q|) = Giikl/Ql) 


(4.6) 
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For a generic term in Q. the transverse momenta k 4 and k B are independent variables, 
while the factor S(k A )d(k B ) turns into l/(k A k B ). Taking moments with respect to s and 
changing integration variables to = ~k A /q A , £ B = k B /q B and A = £,a£,bs/Q 2 gives the 
structure 



d\ N 

T A 


H(\Q 2 ,k A ,k B ) 


(4.7) 


Having evaluated H at some fixed perturbative order, we are interested in the poles of <jn 
at N = 0. There are poles from the small £4 and C, B ends of the integrations over £4 and 
£b. To evaluate < 7 at without losing any powers of 1/1V, we approximate 


J(U -k. 4 ; Qa) - J{ 0, -k a; Q\) = G 1 (k 2 A /Q 2 A ) , 

kg; Qb) —■*■ ^(0, kg; Q B ) — Gi(k B /Q B ) . (4-8) 


The pole associated with the integration over £ A can be thought of as arising from the 
integration over the rapidity of the final state gluon with the largest rapidity. Almost all 
of the momentum fraction £4 is taken by this gluon. Since, in the leading logarithmic 
approximation, this gluon has a rapidity that is much less than that of the quarks in photon 
A, the momentum fraction £ A is negligible compared to 1. Similarly, the pole associated with 
the integration over £g can be thought of as arising from the integration over the rapidity 
of the final state gluon with the most negative rapidity. The poles from the large A end of 
the integration over the function 7i can be thought of as arising from integrations over final 
state gluons with intermediate rapidities. 

The leading logarithmic result can thus be written in the form 

r d 2 k r rl 2 k 

<tn{Qa,Qb) = J Gi(ki/Qi) T N {k\y B ) G.ikl/Ql) , (4.9) 

where .^(k^, kg) is the BFKL function that describes the interaction of gluons with fast 
moving colored systems moving in opposite directions and obeys the BFKL equation U . It 
is normalized so that at order we have 


H 0) (kikD = K 1,5(14_ k |) • 

The solution to the BFKL equation can be written to all orders in a s in the form 

1 ,1/2+ioo dl /k^ 1-1/2 1 


(4.10) 


4v(k 2 ,ky = 


2 ir,/k^k| 4/2-ioo 2-1 lk;J N-a,x{ 7 ) 


(4.11) 


where 


a s = a s C A /vr , C A = 3 


(4.12) 
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The function 7 ( 7 ) is determined by solving the eigenvalue problem for the BFKL kernel and 
is given by 


x(i) = 2 V’(i) - i’h) - V'( i - 7) 


(4.13) 


with being the Euler ^-function. 

The BFKL function ( |4. 1 1|) has poles at N = 0 order by order in perturbation theory. 
Eq. (|4.9|) shows that the poles in the 7 * 7 * cross section are generated from the ones in 
JF/v by integrating the color functions over k±. While these functions are specific to the 
off-shell photon probe, the function T is universal. The same function contributes to the 


small x behavior of the cross sections in hadron-initiated processes [|lfj] via the high energy 
factorization formulas. 

By inserting the representation ( 4.11 ) in Eq. m, and scaling out the dependence on 
the photon virtualities, we get 


a N{Q 2 A, Q 2 b ) — 


2 7r Qa Qb 


/■1/2+10O d ry 
'1/2—ioo 27 Tl 



7 - 1/2 


N - a s x{ 7) 


W(7)^i(l-7) 


(4.14) 


where 1 / 1 ( 7 ) i s defined as the following /c^-transform of the photon color function G\ 

.2 / k 2\7-i 

M 2 ) 


^ 1 ( 7 )=r 

Jo 


d k 2 

~k? 


k 2 


(4.15) 


The explicit expression of the function V\{j) can be determined by using the represen¬ 
tation (|2 .1 7|) and performing the integral transform. The result reads 


Vi ( 7 ) =7raa s 


A (1 + 7) (2-7) F 2 ( 7 ) F 2 (l - 7 ) 

7 (3 — 2 7 ) T(3/2 + 7 ) T(3/2 — 7 ) 


(4.16) 


Eq. (|4.14| ), together with the explicit formulas (fl,13| ), (|4.16|) , gives the leading logarithmic 
result for the moments of the 7 * 7 * total cross section. It sums the 1/N poles to the accuracy 
( a 2 /N ) x (a s /N ) k , for any k. 

The lowest order perturbative contribution, k — 0, can be recovered from the summed 


formula ( |4.14| ) by expanding the denominator (N — a s x(y)) to the zeroth order in a s . 
The simple pole IV -1 is the Mellin transform of unity, and one can check numerically that 
the 7 -integral 


2^Q A Q B ^(Ql,Ql)=r /2+i00 Jl(^ 

Jl/2—ioo Z7TZ \ (cj 


7-1/2 

2 , Vi( 7 ) Vi(l - 7) 

B ) 


(4.17) 


agrees with Eq. ( |2.23|) and Fig. 0. 

In general, multiple pole contributions 1 /N k+X to the cross section are obtained by re¬ 
taining higher orders in the ct s -expansion of Eq. ( |1.14| ). We see that the general structure 
of the coefficients of the leading logarithmic series comes from both 7 ( 7 ) and 1 / 1 ( 7 ): the 
former is a universal function associated with the BFKL pomeron, while the latter describes 
the coupling of the pomeron to a specific physical source. 
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A. Energy dependence 


The total cross section is obtained from Eq. (|4.14| ) by taking the inverse Mcllin-Fourier 
transform ( |4.2| ) . By evaluating the iV-integral from the residue at the pole N = a s x( 7 ), 
one gets 


<r(s, Q \, Q 2 b ) — 


1 /■1/2+ioo d'y 

2 7 r Qa Qb 4i/2 —200 27 ri 



7 - 1/2 



x(l) 


hi( 7 ) Vi(l — 7 ) . (4.18) 


Note that this result depends on two mass scales, the scale / 1 2 in a s and the scale Q 2 in 
the high energy corrections, whose reliable determination would require a next-to-leading 
analysis. We discuss the uncertainties in the leading log result associated with these scales 
in Secs. 0 and [X| 

I 11 the limit s —» 00 the integral (|4.18|) is dominated by the region near 7 = 1/2, where 
the function x has a saddle point. In the saddle point approximation one obtains 


a(s,Q%Q 2 B ) ~ 


|hh(l/2)| S 


5»X( 1/2) 


2-kQaQb ^2 7tx // (1/2) a s In (s/Q 2 ) \Q 9 


(4.19) 


with 


y(l/ 2 ) = 4 In 2 ~ 2.77 , W(l/ 2 ) = 9 7 r 3 a ^ e 2 j /8 , x"(l/ 2 ) = 28 C(3) ~ 33.66 . 

(4.20) 


Eq. (|4. 19|) shows the asymptotic power behavior characteristic of the QCD pomeron, s x 
with A = a s x( 1/2) — 2.65 a s . The pre-factor that determines the normalization of the 
asymptotic cross section, on the other hand, depends on the off-shell photon probe, and is 
controlled by the value of the function Vi ( 7 ) at 7 = 1 / 2 . 

If the two photon virtualities are significantly far apart, corrections of order 
\i\(Q 2 A /Q 2 B )/(a s ln(s/Q 2 )) need to be taken into account when one calculates the large 
s limit. In this case, one finds that the position of the saddle point is shifted to 
7 ~ 1/2 — hi(Q\/Q 2 B )/ (d s x // (l/ 2 ) ln(s/Q 2 )). Apart from corrections in the pre-factor, 
the net effect of evaluating the integral (|4.18|) around the shifted saddle point is to multiply 
the expression in the right hand side of Eq. (|4.19|) by the factor 


exp 


In 2 ( Q\!Q 


2 x"(l/2) ol s hi (s/Q 2 ) 


(4.21) 


The cross section acquires a gaussian modulation in In (Qa/Qb) with a width that grows 

like Ins. 

I 11 the general case, the integral (|4.18|) can be performed numerically. In Fig. |(J we show 
the result as a function of s/Q 2 for a given choice of the values of the photon virtualities and 
the strong coupling. For comparison we also plot the saddle point formula. As the energy 
increases the two curves get closer. However, in the range considered, the sub-asymptotic 
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FIG. 6 . The s/Q 2 dependence of the 7 * 7 * cross section, Eq. (|4.18| ). We take Q 2 A 
a s = 0 . 2 , and divide a by the Born cross section, which eliminates the factor 1 /(QaQb)- 


Q 


2 

B' 


contributions are still significant (about 50% at Q 2 /s ~ 10 _1 , 25% at Q 2 /s ~ 10” 3 ). The 
large size of the corrections to the saddle point approximation can be mainly traced back 
to the fact that the function 1 / 1 ( 7 ) itself is rather sharply peaked around 7 = 1/2. This 
is illustrated in Fig. [7j. where we see that, for instance, for a s = 0.2 and s/Q 2 = 10 2 , the 
width of the pomeron factor is still not small compared to the width of the factor associated 
with the off-shell photon color function. This effect accounts for most of the shift in the 
normalization of the cross section between the asymptotic and exact evaluation of the leading 
logarithmic sum. 


B. Summed results for the asymmetry and the longitudinal cross section 


Using the same method described above for the cross section averaged over the two 
transverse photon polarizations, one can derive summed results for each photon polarization 
(1, 2, L) and for the polarization asymmetry A. Denoting by a, a' and b , b' the polarization 
indices for photons <74 and < 7 b, we write the factorization formula for the moments of the 
polarized cross section 07 / ,bb as 




d 2 k/i r d 2 kf 


7rk 2 4 


?rk| 


(kA, Qa) k B ) G bb (k B , Q b ) 


(4.22) 


The color functions G have been discussed in Sec. fl|. The full solution for the Green’s 
function T of the BFKL pomeron reads [|I|] 
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A 

FIG. 7. The BFKL pomeron factor exp \a s {x{ 7 ) — x(l/2)) ln(s/Q 2 )] and the photon fac¬ 
tor V\ ( 7 ) Vi(l — 7 )/F l 2 ( 1 / 2 ) (normalized to the saddle point) along the contour of integration in 
Eq. (|4.1g| ). We parametrize this contour as 7 = 1/2 + i A. We take a s = 0.2, and show the BFKL 
pomeron factor for two different values s/Q 2 = 10 2 , s/Q 2 = 10 3 . 


^iv(kA, k_e) = 


2 7r \!~k 2 A k| 


E 


m =—00 


P 1/2+200 

d 7 1 

( k2 A 

J 1/2—zoo 

2iri N - a s Xm( 7 ) 

U % 


7 - 1/2 


with 


(4.23) 

Xrn(7) = 2-0(1) -^(7+ |rra|/2) --0(1*7+ M/ 2 ) • ( 4 - 24 ) 

The function T in Eq. (|4.11|) , relevant to the case of the unpolarized cross section, is given 
by the term rri — 0 in Eq. (f4.23| ). 

Consider the product E = a A of the asymmetry times the averaged cross section, with 
the asymmetry A defined in Eq. (|3 .1 1|) . Due to the azimuthal integration only the terms 
with m = +2, —2 in Eq. (|4.23|) contribute. Substituting the color function G 2 (see Eq. (|3T8|) ) 
in Eq. one finds for the N moments 


Z n (Q 2 a ,Q 2 b ) = 


1 (-1/2+ioo rfy / Q 2^' 

16 7 tQaQb J 1/2 —zoo 27 n \Q 2 B/ 


N -a s \ 2 ( 7 ) 


V 2 (j)V 2 (1-j) 


(4.25) 


where V 2 ( 7 ) is defined from the fcj_-transform of the function G 2 analogous to (|4. 1 5|) , and 
has the expression 
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V 2 (j) = 2naa s ( J2 1 


r(i + 7)r(2-7)r( 7 )r(i-7) 

(3 — 2 7) T(3/2 + 7) T(3/2 — 7) 


(4.26) 


By inverse Mcllin-Fourier transformation, one obtains the summed formula for the asymme¬ 
try in the energy space: 


^4(s, Q\, Q 2 b ) — — 


1 


1 


<■t(s,Qa,Qb ) 16 7T QaQb 


X 


/■l/2+i00 ^ ry 

1 1 / 2 -ioo 2m 



7-1/2 



Os X2(7) 


V 2 {l) V 2 (l - 7 ) 


(4.27) 


where cf is given by Eq. (|4.18|) . 

As in the case of the unpolarized cross section, the asymptotic behavior in the limit 
s —> 00 is determined by the saddle point approximation to the integral in Eq. ( 4.27 ). In 


this case we observe a negative power law with the energy s, controlled by the value of \ 2 at 
the saddle point, y 2 (l/2) = —4(1 — In 2) ~ —1.23, indicating that the angular correlations 
between the two photons tend to be washed out when BFKL pomeron exchange dominates. 

The contributions to the photon-photon cross section from longitudinally polarized pho¬ 
tons are given by formulas analogous to Eq. (|4.18| ) in terms of different combinations of the 
functions V\ and V 2 , and the same function y: 


o 


{LL) (s,Q 2 a,Q 2 b ) = 


2 7T Q A Q B 

a (TL \s, Q\, Q 2 b ) = cd LT) (s, Q 2 a , Q 2 b ) 

1 rl/2+ioo d'y (Q 


rl/ 2 +ioo d'y 

1/2—ioo 2m 



7 - 1/2 



0 ‘s x{l) 


V 2 { 1 )V 2 {l-l) 


7 - 1/2 


2 ttQaQ b J 1/2 —200 2ni 



os xh) 


F- 717 ^( 7 )^( 1 - 7 ) 


(4.28) 


Unlike 1 / 1 ( 7 ) i 11 Eq. (|4.16|) , U 2 has simple poles at 7 = 0,1 (instead of double poles), cor¬ 
responding to the non-logarithmic behavior of the color function G 2 at k 2 < Q 2 , k 2 > Q 2 
noted in Sec. 0 . Replacing functions V\ by functions U 2 accounts for the different size of 
the longitudinal cross sections with respect to the purely transverse one. Roughly, one finds 


<7 


« 0.05 a( TT ) , « 0.2 a ^ 


(4.29) 


V. SCALE DEPENDENCE AND UNCERTAINTIES OF THE LEADING 
LOGARITHMIC APPROXIMATION 


The result (|4.18|) for the 7* 7* cross section depends on two mass scales which cannot 
be determined in a leading logarithmic analysis: the mass p 2 at which the running coupling 
a s is evaluated, and the mass Q 2 that provides the scale for the high energy logarithms. 
The former can be thought of as being associated with the integrations over the transverse 
momenta in the loops contributing to higher order diagrams, while the latter stems from the 
longitudinal integrations. A reliable determination of these scales would require a next-to- 
leading order calculation. Lacking this, we provide here qualitative arguments to relate these 
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scales with the physical hard scales of the problem. In Sec. |X] we will use these relations to 
examine numerically the dependence of the cross section on the scale choices. 

A possible choice of the scale /x 2 in the strong coupling is based on the prescription of 
Ref. 0. To apply this prescription, we consider the lowest order gluon exchange contribu¬ 
tion, which we have discussed in Sec. This is given (see Eq. (|2. 1 1|) ) by the integral in dk 2 
of two G factors, each of which is proportional to a s , G = a s (/i 2 ) Q. We first compute the 
quark loop contribution to the gluon propagator, renormalized in the MS scheme: 

n(k 2 ) = A, [ln(k 2 //i 2 ) + C] , (5.1) 


where (3 0 = —1/(6 7r) is the contribution to (3 0 = (33 — 2 Nf) /(12 n) from one quark loop, and 
C = —5/3. Inserting a quark loop into the gluon propagator in the lowest order diagram 
amounts to replacing a s (/i 2 ) by 


a s (n 2 ) 


1 ~(3 0 a s (n 2 ) ln(k 2 e c //i 2 ) 


(5.2) 


(as in m this can be regarded as a contribution to an effective coupling a e //(k 2 )). Now, 
following JlTj we choose the scale /i 2 so that the quark loop contribution vanishes after 
integrating over k 2 : 


f°° dk 2 
Jo (k 2 ) 2 


Stk'VeS) 5(k 2 /<3|) ln(k 2 e c /fi 2 ) = 0 


(5.3) 


This is the same procedure that applies in the case of the abelian theory. Using the repre¬ 
sentation (|2 .1 7|) for G , we are led to calculate an integral of the form 


2 0 c s 


dk" 


An (1 - A A ) k 2 + z A ( 1 - z A ) Q a A B (1 - A B ) k 2 + z B (1 ~ z B ) Q 


In 


k 2 e 1 
v d 2 


(5.4) 


Exploiting the symmetry of the integrand under the transformation k 2 /Qa Qb —* QaQb/ k 2 
and interchange of the variables zn X B , z B <->■ An, one can show that the condition ( p.3| ) 
is satisfied by 


d 2 = c i1 QaQ 


B 


c „ = e-'/ 3 


(5.5) 


Note that, when the two photon virtualities Q\ and Q 2 B are far apart from each other, the 
prescription 0 picks out a scale for a s which is neither of the order of the big virtuality nor 
of the order of the small one, but rather is proportional to the geometric mean \JQ 2 a Q 2 b . The 
value of the proportionality coefficient in Eq. ( |T5|) is specific to the subtraction scheme 
(MS) chosen to define the quark loop insertion. 

The argument given above applies to the factors of a s that appear in the Born ap¬ 
proximation to the high energy cross section. The result (|4.18|) , however, also contains a 
dependence on the running coupling through the higher order factor s“ sX associated with 
the solution of the BFKL equation. For the scale in a s in this case one does not have such 
a simple argument as the one described above. For the numerical estimates in Sec. |X| we 
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will make the assumption that the same value of p 2 also controls the running coupling in 
the BFKL factor. In Sec. |X| we will check the numerical effect of varying this scale. 

We now consider the mass Q 2 that provides the scale for the large energy logarithms 
(a s In (s/Q 2 )) k (see Eq. (|4.1| )). To estimate this scale, we observe that the rapidity of gluons 
exchanged in the rungs of the BFKL ladders should lie between the rapidity i/a of the quark 
Pa (produced by the photon g^) and the rapidity ys of the quark ps (produced by the 
photon g#). This gives rise to integrations over the rapidity intervals 


rvA 

/ dy = y a-Vb 

J VB 


(5.6) 


The logarithms of the energy s are generated precisely by these integrals. Estimating the 
size of the rapidity intervals thus allows us to estimate the scale Q 2 . 

Expressing the “+” and ” momentum components of the quark pA as 


Pa = z a qj 


Pa = 


Pj LA 

2 pi 


(5.7) 


we write its rapidity as 


SA = i.„ (rtWihd 2 ^)- 

2 \p- A ) 2 


r» 2 

P_L A 


(5.8) 


Similarly, for the quark ps we have 


-HS, 


=-ii,d 24( « 


2 '~b) 2 


pI b 


(5.9) 


Taking the difference between these rapidities gives 

Z A Z B S 


VA~yB = hi 


\P±a\ |P±b| 


(5.10) 


The average transverse momenta p_la, Pub carried by the quarks are of the order of the 
photon virtualities Qa, Qb■ For the longitudinal momentum fractions za, zb, we assume 
a typical maximum value of the order z max ~ 0.1 in the high energy region. Using these 
estimates, we obtain 


Va-Vb = In 


( g 4ax \ 
\Qa Qb J 


(5.11) 


We thus identify the scale 

Q 2 — C Q Qa Qb 


CQ = V4ax = !0 2 


(5.12) 
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VI. LIMITATIONS ON THE BFKL POMERON APPROACH AT VERY LARGE 

ENERGY 


The transverse-momentum integrations in the factorization formula ( |4.9| ) are dominated 
by values of k^, of the order of the photon virtualities. (For the purpose of this section 
we will take the photon virtualities to be equal, Q\ = Q 2 B .) As a result, for sufficiently off 
shell photons the dominant contribution to the cross section comes from short distances, 
and the evaluation of Eq. ( [4.9| ) gives rise to a finite result in perturbation theory, ffowever, 
there are limitations on the perturbative treatment that are intrinsic to the BFKL equation. 
These limitations come from the region of very high s. Although an accurate understanding 
of these effects is an open problem ||T 8 [ that goes beyond the scope of this work, one can 
nevertheless make some rough estimates. We discuss them in this section. In order to keep 
the notation simple, we do not distinguish here between the scales QaQb, Q 2 = cqQaQb, 
and /i 2 = c^QaQb, calling all of these simply Q 2 . 

It is known from the structure of the BFKL equation that, even if the incoming k 2 4 , k^ 
are large, say, k.\, k^ ~ Q 2 , the typical transverse momenta in the gluon ladders contributing 


to the function jF^k^, k s ) (see Eq. ([4.11|) ) may diffuse away from Q 2 as the energy becomes 
very large. The diffusion coefficient in In k 2 can be read directly from the asymptotic solution 
to the BFKL equation (or, equivalently, from the exponential term (j4.21| ) in the 7 * 7 * cross 
section) and is given by 


Pdiff = 2x"(l/2)a s ln(s/Q 2 ) . (6.1) 

As a result, the distribution of the transverse momenta in the BFKL ladders is a gaussian 
in Ink 2 centered around Q 2 with a width proportional to pdiff [0- With increasing s the 
distribution broadens, and one becomes sensitive to the region of small transverse momenta. 

A self-consistency check of the perturbative treatment requires that the contribution 
from transverse momenta of the order of A.qcd be suppressed. In order to keep away from 
momenta of this order one has to have 


In 2 (<? 2 /A 2 )/ (2 x"(l/2) a. In( s /« 2 )) > 1 , (6,2) 

Identifying a s with the strong coupling evaluated at the scale Q 2 , a s ~ (/3 0 ln(Q 2 /A 2 )) 1 , 
one gets 

a s (Q 2 ) In {s/Q 2 ) < d/a 2 (Q 2 ) , d = tt/ (2 y"(l/ 2 ) /3 2 C A ) « 1/30 . (6.3) 

For any given Q 2 , this can be read as an upper bound on the domain of energies s in which 
we expect the perturbative approach based on the BFKL pomeron to be reliable. Observe 
that, for small values of a s , this is not so stringent a constraint. However, because the 
BFKL function y has a large second derivative at the saddle point, the numerical value of 
the coefficient <7 is small. Therefore, the limit ( |6.3| ) may become relevant unless Q 2 is very 
large. 

The BFKL equation is also known to give rise to violation of the unitarity bound at 
asymptotically large energies. The growth of the cross section predicted by the BFKL 
equation cannot continue indefinitely, and unitarity corrections must arise to slow it down. 
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Roughly, these effects are expected to become important when the calculated cross section 
is bigger than the naive geometrical cross section 1 /Q 2 . A careful discussion may be found 
in Ref. [:2(J. The simplest estimate, a < 1 /Q 2 , yields, using Eq. (|4. 19|) , 


<X 


2 7T x' / (l/2) a s in (s/ (Q 2 )) 



«3X(l/2) 


<1 


(6.4) 


where we have neglected factors from the term in TA (1/2) in Eq. (|4.19|) . Assuming, as we 
did before, a s to be evaluated at the scale Q 2 , we can write 


a s (Q 2 ) In {s/Q 2 ) < c 2 ln(l /a 2 {Q 2 )) + A 


c 2 = t r/ {C A x{ 1 /2)) ~ 1/3 


(6.5) 


In the term A we collect contributions arising from the factor in the square root in Eq. (El 
as well as from the factor V\ in Eq. (|4.19| ). 

The different functional dependence on a s in the right hand sides of the inequalities 
o and ( |6.5| ) implies that, for small enough values of a s , the unitarity limit ( |6.5| ) is more 
stringent than the diffusion limit (|6.3|) . This suggests that for sufficiently high Q 2 it should 
be possible to study unitarization in a purely perturbative context [[2D] . 

On the other hand, the coefficients c± and c 2 are significantly different in size. This 
may make the two bounds (0) and ( |6.5| ) rather comparable for moderate values of Q 2 . 
In addition, inspection of Eqs. (|4.19|) , ([O]) suggests that the term A in Eq. ( |6.5| ) is not 
necessarily negligible, and may contribute to push the onset of unitarity corrections further 
away. 

As to the impact on experimental studies at future e + e~ colliders, we observe that 
unitarity corrections should set in when the cross section has grown to be much larger than 
the Born value ~ a 2 /Q 2 . At a future e + e~ collider one may thus expect to see the rise of 
the cross section with s/Q 2 . Possibly, one may see this rise slow as unitarity corrections 
become important [10] . 


VII. VIRTUALITY DEPENDENCE AND RELATIONSHIP WITH DEEPLY 

INELASTIC SCATTERING 


If we let Q\ be much larger than Q\ in Eq. (|4.18|), we obtain a result that describes 


small x deeply inelastic scattering from a transversely polarized photon [2T| whose virtuality 
Q 2 b is sufficiently large to allow the use of perturbation theory to analyze its decomposition 
into quarks. The structure function Fi{x, Q 2 ) (which is the same as F T ) is related to the 
virtual photon-photon cross section by 


Fi(x,Q 2 a ) = 


(2n) 2 a 


Q a Q Ai Q 


with 


s = 


Ql 


X 


B) 


(7.1) 


(7.2) 


Thus Eq. ([4.18|) gives a leading logarithmic summation for F x at small x: 
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(7.3) 


Fi(x,Q 2 a ) 


1 

(27r) 3 a 


/•1/2+zoo 

' 1 / 2 —ioo 


d 7 
2 m 


oiy (Qa_Y‘ xM 
QbJ {zQ 2 J 


The scale Q 2 is not fixed by a leading logarithmic calculation. A sensible choice for Q 2 
based on a qualitative argument has been discussed in Sec. 0. For the purpose of the 
present section we simply leave it undetermined. 

Consider the perturbative expansion of Eq. (17.31 ), 


Fi(x,Q 2 a ) 


1 

(27t ) 3 a 


/“1/2+200 

'1/2—200 


d 7 
27rz 



1 + a s x{l) In 




Ei( 7 )Ei(1- 7 )- 


(7.4) 


The functions 1 / 1 ( 7 ), V\(l— 7 ) and 7 ( 7 ) have poles at integer values of 7 . Thus this expression 
contains contributions proportional to (Q\)° times logarithms, (Q ^) _1 times logarithms, 
and so forth. We extract the leading twist contribution, the part proportional to (Qa)° 
times logarithms, by rewriting the integral as an integral over a contour C that encircles the 
singularity at 7 = 0 plus an integral over a contour from —1/2 — zoo to —1/2 + zoo. The 
integral over the contour C is the leading twist contribution. We discard the integration over 
the contour from — 1/2 — zoo to — 1/2 + zoo, which contains the higher twist contributions. 
Thus the leading twist contribution to iq is 


FiMa) = 


1 [ dj_ _ 

(27r) 3 a Jc 2m 


Res 



1 + x(l) In 


Q 


x(P 


(27r) 3 a 7=0 \Q 2 b ) 



1 + a 8 x{l) In 


'Q\_ 

xQ 2 


+ ••• Vi( 7 ) V\{1 — 7 ) 

+ ... Ri( 7 )Ri( 1 - 7 ). (7.5) 


It is straightforward to determine the perturbative coefficients. The function V\(j) has 
a double pole at 7 = 0, as does the function Vi (1 — 7 ). Thus there are three powers of 
In (Q\/Q 2 b) a t leading order in a s , which is a 2 since V\ oc a s . We find 


f l.t. 


(x,Qa) 


4 a a 2 
9 7T 3 




+ - In 2 
3 




+ 



f C(2)) +o(«,) 


. C(2) ^ 1.64 . 



(7.6) 


If we restore the summation of the leading logs of 1/x, we obtain 


Fi(*,Qa) = 


1 /• h (Q\ y (q 

( 2n) 3 a Jc 2m \Q%) 


x(7) 

El( 7 ) Vi(l — 7 ) 


(7.7) 


This is not a simple function. However, one can easily perform the integration numerically. 

One interesting result is that we observe what used to be called “precocious Bjorken 
scaling”. That is, the leading twist approximation begins to be accurate at values of Q\ 
that are not really very large compared to the only other scale in the problem, Q\. We 
illustrate this in Fig. working at order a 2 . We compare the leading twist approximation 
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that applies for Q\ Q 2 Bl the leading twist approximation that applies for Q 2 B Q\, and 
the full result. First we define xa = Q\/ s and plot 

<tl.t.-a(s, Q 2 a,QI) = a F i L ' T ' ( x -4, Q\) (7-8) 

versus Qa/Qb • Then we define xb = Q' B / S and plot 

o-l.t.—b(s, Q 2 Ai Q\ ) = (2 ^ 2 a Ft T ix B , Q%) (7.9) 

Vb 

Finally, we plot cr(s, Q% Q 2 B ) at the leading order in a s without approximation. We see that 
the leading twist approximation for Q\ Q B is quite good down to Q\ = Q 2 B , while the 
leading twist approximation that applies for Q 2 B Q\ is quite good down to Q 2 B = Q\. In 
fact, at Q\ = Q 2 b , both approximations are quite good. A similar behavior is observed if we 
include higher orders in a s . 



Q a / Qb 

FIG. 8. The full and the leading-twist approximation to the virtual photon cross section. We 
plot the cross section divided by 16a 2 a 2 (Y, q e 2 ) 2 /(QaQb), as in Fig. [|. 


VIII. REGGE FACTORIZATION 


The QCD result for 7* 7* scattering can be compared with expectations for the structure 
of the high energy cross section based on traditional Regge theory . In Regge theory, to 


analyze the elastic scattering of particles A and B one considers the singularity structure of 
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the amplitude A(s,t) in the complex angular momentum plane. The simplest case is a pole 
in the angular momentum j plane located at a position dependent on t, 1 / (j — a(t)). One 
then obtains the asymptotic behavior s a(t > for s —> oo and t fixed, and the amplitude takes 
the factorized form 


2l(s, t) ~ P A (t) s“ (,) P B (t) 


(8.1) 


Here (5a, (5b are functions of the transferred momentum, associated respectively with the 
couplings of particles A and B to the reggeon whose trajectory is a{t). 

For the case of the 7* 7* total cross section, this would correspond to the structure 


a. 


7 7 


(s, Q\,Ql) ~ Pa( f = 0; Q 2 a ) P B (t = 0; Q 


B) 


( 8 . 2 ) 


where we have used the optical theorem to relate the total cross section to the imaginary 
part of the forward elastic amplitude, and we have let ( 5 a and ( 5 b depend on the photon 
virtualities in the case of off-shell photons. 

However, it has been long known, from various phenomenological and theoretical consid¬ 
erations |E|, that this structure cannot be exactly true, and strong-interaction scattering 
at high energy has to have a more complicated singularity structure than a pole, such as 
moving or fixed cuts. In this case, one does not expect the factorized form for the cross 
section to hold. 

If we now turn to the QCD result (see Eq. ( f 4 . 14 | )), we may ask what kind of singularities 
the amplitude has in the angular momentum plane. The BFKL pomeron is known to give 
rise to a (fixed) branch point singularity [|l|]. To see this, let us consider the moments of 
the 7* 7* cross section. At leading level, we are allowed to identify the moment N with 
the complex angular momentum j 0. Eq. ( [ 4 . 14 ] ) is written in terms of an integral in 
the complex 7-plane. The integrand has the pole 1 / (N — 07 7(7)) and one factor of V for 
the coupling of the gluon system to the quarks in each photon. Thus the integrand has 
a factorized structure. To understand the angular-momentum singularity structure, one 
should see what becomes of the pole N = 077(7) after 7-integration. 

For Q\ > Q B , the 7-integral is well approximated by the residue at the rightmost pole 7 
to the left of the integration contour. As we have seen in the previous section, numerically 
this approximation turns out to be fairly good down to values of Q\/Q 2 B just above 1. In 
this approximation one gets 


°n{Q 2 a, Q 2 b) 


2 ttQ 2 


C N exp 7 at In {Q\/Q 


B) 


C N = 


Vi(7 n) Vj(l — 7 n) 
UsX'ilN) 


( 8 . 3 ) 


The leading pole 7 n is determined by the equation N — 077(7) = 0. At A^ = 077(1/2) it 
has a square-root branch point singularity: 


7jv ~ 2 - 


2 (.N - N 0 ) 


N —> IVo = 07 7 ( 1 / 2 ) 


\| <XsX"{ 1 / 2 ) ’ 

Correspondingly, the residue Cn is also singular: 

C N ~ |U(l/2)| 2 fo2a s x"(l/2) (N-N a ) , N^N„ . 


(8.4) 


( 8 . 5 ) 
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Therefore, the QCD result implies a branch point rather than a simple pole in the 
angular momentum plane. As a consequence, we do not obtain a Regge-factorized form for 
a(s,Q 2 A ,Q%). 

It is of interest, however, to see by how much this factorization is violated. We first 
consider the asymptotic formula ([4.19|) , obtained by using the saddle point approximation 
around 7 = 1/2. In this case, provided the scale Q 2 is fixed by Q 2 = cq QaQb, as discussed 
in Sec. |V|. and the QCD coupling a s is fixed, an approximate form of Regge factorization 
is recovered. The piece which violates this factorization in Eq. (|4.19|) is proportional to the 
square root of a logarithm, and is therefore a slowly varying function. A more substantial 
source of factorization breaking comes in when one takes into account the correction due to 
the position of the saddle point drifting away from 1/2, see Eq. (|4.21|) . 

For the exact leading log result ( |4.18|) , a possible way to quantify the deviation from the 
Regge-factorized behavior is to look at the ratio 


R = 


a. 


7 7 


(s,Q 2 a,Q 2 b)Y 




a. 


( 8 . 6 ) 


7 7 


If Regge factorization holds, this quantity should be equal to 1. On the other hand, from 
Eq. ( |T| ) we see that R goes like R ~ Q 2 B /Q 2 A for Q 2 A Q%, and R ~ Q\/Q\ f° r 
Q 2 b Q\. That is, for Qa and Qb sufficiently far apart the Regge-factorized form breaks 
down. 
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FIG. 9. The Qa/Qb dependence of the ratio R defined in Eq. (jOi|), at different values of s/Q 2 . 

In Fig. ^j. we plot R as a function of the ratio of the photon virtualities Qa/Qb f° r 
different values of the energy. It is interesting to observe that, for typical parameter values, 
















R varies by not more than 40% when Qa/Qb varies from 0.1 to 10 . We take this as an 
indication that an approximate Regge factorization holds numerically for the exact integral 
( |4.18| ) if Qa/Qb is not too large or too small. 


IX. SOFT SCATTERING AND HARD SCATTERING 

The calculation of the high energy photon-photon cross section discussed so far is a 
perturbative calculation, based on the dominance of short distances for large photon virtu¬ 
alities. When the photon virtualities decrease, one goes out of the region of validity of the 
perturbative approach. As the photons go near the mass shell, the high energy scattering 
process is expected to become dominated by soft interactions. In this regime one is not able 
to calculate in QCD, and in order to have a (phenomenological) description of the cross 
section, one rather has to rely on models for strong-interaction scattering based on Regge 
theory. 

For on-shell photons, the Regge factorization hypothesis allows us to relate the photon- 
photon total cross section to the photon-proton and proton-proton cross sections, as follows 

(Tryry ~ O'ryp Oryrp / (7p p . (9T) 

Assuming the values o lp rs 0.1 mb, o pp ~ 40 mb, one gets a 77 ~ 250 nb. For virtual 
photons with small Qa and Qb, the fall-off of the cross section can be estimated from vector 
meson dominance: 


o. 




( K V 

( M ' \ 

\M 2 p +Q 2 a) 

\M/ + Q 2 b J 


a, 


77 


(9,2) 


As the photon virtualities increase, the cross section, instead of continuing to fall like 1 /Q 8 , 
should begin to fall more slowly. At large photon virtualities (of the order of a few GeV, 
or bigger), it should go over to the perturbative scaling behavior in Eq. (|4.18| ), cr oc 1/Q 2 at 
fixed s/Q 2 . 

Note that the 1/Q 2 behavior could not be obtained in the framework of the Regge 
factorization (|9.1| ) even if one assumed perturbative scaling in each one of the photon cross 
section factors, that is, even if one assumed o lp oc 1 /Q 2 . This is the counterpart (at the 
level of hadronic cross sections) of the effect of the deviation from unity observed for the 
ratio R in the previous section (see Eq. ( | 8 . 6 |) ). In fact, experimental data on the yp cross 
section for large photon virtuality are now available in the region of high energies from 
the measurements of small-x deeply inelastic scattering at HERA. The above observation 
amounts to saying that, even if one used the data for cr 7P , the relation (|9.1| ) would not lead 
to the correct perturbative QCD result for the 7 * 7 * cross section. 

I 11 Fig. [[(] we show a log-log plot of the curves corresponding to the soft and perturbative 
formulas for the Q 2 -behavior of the cross section [|j] . For the former, we use Eqs. (0-(ED> 
and for the latter we take the Born approximation to Eq. (|4.18|) . I 11 this plot we are interested 
in emphasizing the dependence of the cross section on Q 2 at fixed s/Q 2 . For this reason we 
limit ourselves to the lowest order formulas, and do not include the higher order summation 
of the ln(s/Q 2 ) terms that give an enhancement at large energies. It is understood that 


29 














FIG. 10. 
order, with 


Q 2 -behavior of the vector meson dominance and perturbative cross sections in lowest 

Qa = Qb = q 2 - 


such high-energy corrections affect both the soft and perturbative curves, giving rise to 
“soft-pomeron” and “hard-pomeron” effects in the two cases. 

The region of intermediate values of Q 2 in Fig. |10] (Q of the order of 1 GeV) is where the 
transition from the soft-scattering regime to the hard-scattering regime is expected to occur. 
The mechanism through which this happens is not theoretically under control at present, 
and one may consider trying to estimate the cross section in this region by interpolating 
between the two curves. 

In the next section we will study the prospects for investigating high energy 7 * 7 * scat¬ 
tering at e + e~ colliders, and we will discuss which regions in Q 2 in Fig. 10 can likely be 
accessed experimentally at LEP200 and a future e + e~ collider. 


X. NUMERICAL RESULTS FOR THE ELECTRON-POSITRON CROSS 

SECTION 

The cross section for high energy virtual photon scattering can be measured in e + e~ 
collisions in which the outgoing leptons are tagged. The cross section for the electron- 
positron scattering process is obtained by folding the 7 * 7 * cross section with the flux of 
photons from each lepton. Consider the four-fold differential e + e~ cross section averaged 
over the angle between the lepton scattering planes, Eq. (|1.2|) . To get an estimate of the 
rates available to study BFKL effects in virtual photon scattering at e + e~ colliders of the 
present and next generation, we integrate this cross section over a region 1Z determined by 
cuts that we discuss below: 
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( 10 . 1 ) 


a = 


dxAdx B 


dQ\ dQ 2 B 

Q\ Ql 


Qa Q% da ( ' e+e ] 
dxA dx B dQ\ dQ 2 B 


We choose 

i) Qa > Q min, Qb > Qmm, where Q min is a few GeV, in order that the coupling a s be 
small, and that the process be dominated by the perturbative contribution; 

ii) Xa x B s ee > bQaQbi in order that the high energy approximation be valid. We 
discuss the parameter k below. 

Note that, with these criteria, the photon virtualities Q 2 lie in a range Q 2 min < Q 2 <C s in 
which the equivalent photon approximation (Eq. ( |1.2|) ) is expected to work fairly well. On 
one hand, kinematical corrections of order Q 2 /s are suppressed in this range. On the other 
hand, contributions of order m 2 /Q 2 to the e —> eq splitting process can be neglected. 

To choose a value for the parameter k, we compare the gluon exchange contribution 
with contributions that are suppressed by a power of s H$|]. We consider first the two gluon 
exchange graph, for which a ~ s° for large s. Taking the case Qa = Qb = Q, we have from 
Fig. 1 


cr 


( 0 ), 


s,Q 2 


0.9 x 16 


a 2 a 2 




e 2 
q z-q 


Q 2 


(l + 0(Q 2 / s)) 


( 10 . 2 ) 


Next, we consider the leading order (electromagnetic) contribution to 7 * 7 * —*► qq, occur¬ 
ring via quark exchange, for which ~ 1/s. In Appendix [B] we report results for this 
subprocess. For Qa = Qb = Q and large enough s, the corresponding cross section is well 
approximated by the formula: 

a ( 9 ) (s,Q 2 ) w -— ^ ^ q q (ln(s/Q 2 ) - l) (l + 0(Q 2 /s)^j . (10.3) 

Demanding that the gluon exchange graph give a larger contribution than the quark ex¬ 
change graph leads to the requirement 


g 1 > 1 8nE q e 4 q ^ 

Q 2 (hi(s/Q 2 ) - 1 ) ~ a 2 o g x 16 (^e ^) 2 a s 


(10.4) 


For typical perturbative values, a s ~ 0.2, we get s/Q 2 > 10 2 . We will therefore use 10 2 as 
a standard value for k. 

We note that for s/Q 2 > 10 2 , one is surely entitled to drop terms in the gluon exchange 
graphs that are suppressed by powers of Q 2 /s, as we have done. 

We thus compute the integrated rate a in Eq. (| 10 . 1 |) using the results given in Sec. |IV] for 
the photon-photon cross section, and setting the scales in the running coupling and in the 
high energy logarithms according to the prescriptions discussed in Sec. 0. The dependence 
of cr on the lower bound Q ^ in on the photon virtualities is shown by the “summed” curve 
in Fig. pT] for the energy of a future e + e _ collider. Fig. [12] shows the cross section for the 
LEP collider at CERN operating at a/s = 200 GeV. 

The dashed and solid lines in Figs. [11] and [12] correspond to the result of using, respec¬ 
tively, the Born and the summed expressions for the photon-photon cross section. At the 
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FIG. 11. The Qmin dependence of the integrated rate a, Eq. ( 10.1 ), for s/s = 500 GeV. We 
take k = 10 2 . We set the scales /r 2 and Q 2 according to the prescriptions given in Eqs. ( |T5| ) and 
( 5.12 ). The solid curve represents the full leading log summation, while the dot-dashed curve shows 
the Born result. The dotted curve shows the contribution to the fully summed result coming from 
transversely polarized photons. 


values of V® considered in the figures, summation effects enhance the rates significantly in 
the range of Q min of a few GeV. As Q min increases, lowest order perturbation theory gets 
closer and closer to the fully summed prediction, as a result of both a s becoming small and 
the phase space closing up for the high energy logarithms. 


In Figs. |TT] and [12], we also plot separately the contribution to the cross section from 
purely transverse photons, that is, the contribution from the term in in Eq. ( |1.2| ). We 
see that this contribution accounts for about three quarters of the full cross section. 


For values of the cuts Q m \n = 2 GeV, n 

= 10 2 , we find 


a ~ 1.5 pb 

(Vs = 200 GeV) 

(10.5) 

at LEP200 energies, and 



cr ~ 12 pb 

(Vi = 500 GeV) 

(10.6) 


at the energy of a future collider. These cross sections would give rise to about 750 events 
at LEP200 for a value of the luminosity L = 500 pb -1 , and about 6 x 10 5 events at a/s = 
500 GeV for L = 50 fb" 1 . 

The choice of the cuts that can realistically be implemented is affected by experimental 
constraints. In particular, the lowest photon virtualities that can be reached are limited 
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FIG. 12. Same as in Fig. |Tl] for y/s = 200 GeV. 


by the angular acceptance of the detector, according to the relation Q 2 « (1 — x)E 2 9 2 , 
where E is the beam energy, 9 the angle of the tagged lepton, and x is the momentum 
fraction of the emitted photon. In this situation, the value Qmm = 2 GeV, for which the 
rates (|10.5|) , ( |10.6| ) are given, implies detecting leptons scattered through angles down to 
about 20 mrad at LEP200, which is close to the range of the current luminosity monitors 
at the LEP experiments |^3|. For a future 500 GeV collider, Q m i n = 2 GeV corresponds to 
a minimum angle of about 8 mrad. It appears that working down to such an angle will be 
difficult but not impossible [^4|. If instead we take Q min = 6 GeV, the minimum angle is 
24 mrad. Then the cross section is about 2 x 1CV 2 pb, corresponding to about 10 3 events. 

As stated earlier, the numerical results given above depend on the choice of the scales in 
a s and in the high energy logarithms that enter the photon-photon cross section. For the 
calculations described above, we have used the prescriptions given in Sec. 0. Different scale 
choices are possible, and they would affect the predictions at the next-to-leading logarithmic 
order, which is beyond the present theoretical accuracy. We can use the variation of the 
results with the scale choices to get an estimate of the uncertainties associated with unknown 
sub-leading corrections. We can vary the two scales p 2 and Q 2 (see Eqs. (|5.5|), (|5.12|) ) 
independently. An illustration of this is reported in Fig. 0- Here we compare the result of 
Fig. |TT] with the curve obtained by multiplying the scale in a s by a factor of 4, /i 2 —> 4 x /i 2 , 
and the curve obtained by reducing the scale in the high energy logarithms by a factor of 
4, Q 2 —> Q 2 /4. The band between these two curves indicates that the uncertainty on the 
leading logarithmic result is fairly large, and emphasizes the need for improving the accuracy 
of the calculations at high energy. 

The plots of the cross section versus Q 2 lin shown in Figs. 0 and [12] illustrate the expected 
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FIG. 13. Estimate of the uncertainty on the leading logarithmic result for the rate a. The solid 
curve is the summed result shown in Fig. 11. The dot-dashed curves summarize the variation of 
this prediction as a result of varying the scales in the strong coupling (y? —> 4/r 2 ) and the high 
energy logarithms (Q 2 —> Q 2 / 4). 


dependence of the photon-photon cross section on the photon virtualities. If we fix Q n 
can look at the dependence on the photon-photon c.m. energy \/I. It is useful to use 


we 


S = Xa'XbS 


and the photon-photon rapidity 



as variables instead of x A and xg. Then we define 
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Using Eq. (S3). we can write dcr/(dins dy) as 
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where we omit three similar terms. We see that da/ (d In sdy) is very directly related to the 
77 cross section. 

We plot da/(d In s dy) at y = 0 in Fig. [L4] for y/s = 200 GeV and in Fig. |1^ for y/s = 
500 GeV. Here we choose Q 2 lin = 10 GeV 2 and k = 10 2 . In each case, we show a curve 
for the Born level cross section and another for the full BFKL cross section. We also show 
the cross section arising from the scattering of (transversely polarized) photons via quark 
exchange instead of gluon exchange. We see that, with our choice of cuts, quark exchange 
scattering is suppressed. 




30 100 500 

Vs [GeV] 


FIG. 14. The cross section da/(dln s dy), Eq. (|I0.9|) , at y = 0 for y/s = 200 GeV. The solid 
curve is the summed BFKL result. The dot-dashed curve is the Born result. The dashed curve 
shows the cross section arising from the scattering of (transversely polarized) photons via quark 
exchange. The cuts are Q 2 in = 10 GeV 2 and k = 10 2 . 


For a /§ < 50 GeV the cross section shows a strong dependence on the cut hQaQb < s. 
With our choice of k — 100 and with Qa,Qb > Qmin = \/l 0 GeV, the cross section is 
forced to vanish for \/§ <31.6 GeV. As \/I increases, the effect of this cut on the Qa and 
Qb integrations becomes less and less important. Since a 1 * 1 *(s,Q\,Q%) is independent of 
s at Born level, the cross section begins to flatten out as \/I increases to about 100 GeV. 
For larger values of \/I, the Born cross section decreases because of the influence of the 
photon flux factor xaP(xa) xbP(xb) = (s/ s)P([s/ s ] 1 / 2 ) 2 . For the summed BFKL curve, 
the growth of a^*^*(s,Q\,Q%) overcomes the effect of the photon flux factor, so that the 
cross section rises with y/§. 

The curves for y/s = 200 GeV and y/s = 500 GeV are similar. The main difference is 
that at y/s = 500 GeV there is more available range for y/l. 

Our plots are for y — 0. The available range of y is \y\ < \ ln(s/s). Thus the cross 
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FIG. 15. Same as in Fig. |lj but for sjs = 500 GeV. 


section integrated over y goes to zero as —»► y r s. 

We see from the results presented above that at a future e + e _ collider it should be possible 
to probe the effects of pomeron exchange in a range of Q 2 where summed perturbation 
theory applies. One should be able to investigate this region in detail by varying Qa , Qb 
and s = x a % B s ee independently. At LEP200 such studies appear to be more problematic 
mainly because of limitations in luminosity. Even with a modest luminosity, however, one 
can access the region of relatively low Q 2 in the graph of Fig. [HI] if one can get down to small 
enough angles. This would allow one to examine experimentally the transition between soft 
and hard scattering. 

We now move on to the angular distribution for the e + e~ scattering cross section, and 
consider the asymmetries Ay, A 2 introduced in Eq. (0)- As pointed out in Sec. 0. Ai is 
zero at leading order. On the other hand, A 2 is given by an equivalent-photon formula in 
terms of the asymmetry A for the 7 * 7 * scattering process discussed in Sec. [IV B| . This reads 


Q\Q\dA e+e } 
dxA dxs dQ\ dQ 2 B 

( q/ \ 2 l _ % A 1 B 

— j A^i* 'y* (xa % B S 5 Q A"! Q b)^ 7* 7*(*^ A % B S 1 Q A"! Q i?) 5 (10.11) 


with A and a being given in Eqs. (|4.27|) , (|4.18|) . 

We can use the same cuts discussed earlier in this section to integrate Eq. ( |10.11|) , and 
thus define 


A 2 


I f a a d ^ A d ®B a Q 2 AQ 2 B d (y {e+e ] 
a Jn A B Qa Q 2 b 2 dxA dx B dQ\ dQ 2 B 


( 10 . 12 ) 
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where a is the integrated rate in Eq. (|10.1|) . Performing the integral numerically, we find that 
the asymmetry A 2 is very small. As noted in Sec. |TV| . the role of the summed BFKL terms 
is that of reducing the magnitude of the asymmetry with respect to the Born order result. 
At a 500 GeV collider, in the range of the angular and energy cuts previously described, we 
find A 2 ~ 10~ 3 . We observe that spin effects in photon-photon scattering at high energies 
are interesting, but the predicted asymmetries are either zero or small. 


XI. CONCLUSIONS 

Understanding the behavior of high energy hadron reactions from a fundamental perspec¬ 
tive within QCD is an important goal of particle physics. As we have shown in this paper, 
virtual photon scattering 7 *(Q\) +7 *(Qb) hadrons at high energies, s » Q\,Q B , pro¬ 
vides a remarkable window into pomeron physics. The total cross section can be studied as 
a function of the space-like mass of each incident projectile. Most importantly, the process 
can be investigated in the regime where the photons both have large virtuality, so that one 
can use the framework of perturbative QCD. 

Compared to tests of the QCD pomeron behavior based on deeply inelastic structure 
functions, the measurement of the total cross section for sufficiently off-shell photons is free 
from the long-distance ambiguities related to the structure of the hadronic target. On the 
other hand, unlike tests based on associated jet production in lepton-hadron or hadron- 
hadron collisions, the 7*7* measurement is fully inclusive and therefore it does not depend 
on specifying the details of the final state. 

The scattering of highly virtual photons can be described as the interaction of two 
incident color singlet qq pairs of small transverse size interacting through multiple gluon 
exchange. We have studied this reaction both in the Born approximation (corresponding to 
two-gluon exchange) and also with the inclusion of the higher-order summation encompassed 
by the BFKL equation. The cross section at high energies and large virtuality takes a 
factorized form in transverse coordinates. However, it does not factorize simply into separate 
functions of Q 2 A and Q B , which reflects the cut structure of the BFKL pomeron in the 
complex angular momentum plane. We have also examined the background contribution 
from quark exchange, a process which is power-law suppressed at high energy. 

According to this analysis, the 7*7* cross section falls off at high virtuality only as 
1/Q 2 , where Q 2 ~ max{Q\, Q B }. The rate for sensitive tagged-lepton experiments at high 
energy e ± e~ or / n ± /r _ colliders is thus not negligible. In particular it appears that the 
main features of the perturbative QCD predictions, such as the energy dependence, the 
factorization properties of the cross section, the scaling laws in Q\, Q 2 B , as well as the 
polarization and azimuthal correlations can be tested in detail at a high-energy and high- 
luminosity next linear collider. We have also found that an interesting first look at virtual 
photon scattering can be obtained from the tagged lepton events measured in the luminosity 
monitors of present experiments at LEP200. 

More precisely, we estimate that, in the region of photon virtualities where summed 
perturbation theory is expected to apply, there should be several hundred events at LEP200, 
and about 10 5 events at a future 500 GeV collider with an integrated luminosity of 50 fb _1 . 
We also find that the enhancement due to BFKL pomeron terms over the Born cross section 
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is sizable, and should be visible particularly in the s-distribution of the cross section, with 

S Xa %B ^ee* 

The dependence of the cross section on the photon virtualities Qa and Qb is perturbative, 
and can be predicted in the framework of the BFKL equation. These predictions can be 
tested by measuring the angles of the recoil leptons. Both the case in which the two photon 
virtualities are varied together (Qa ~ Qb) and the case in which they are kept far apart 
(Qa ^ Qb) are of interest. In the second case one gets to observe the structure function of 
a virtual photon at small Bjorken-x. 

The spin structure is rich, but hard to observe. Most of the observable cross section 
comes from the scattering of two transversely polarized photons. For this part of the cross 
section, there is an asymmetry in the angular distribution of the outgoing leptons, but this 
asymmetry is less than 1%. 

In the region of low photon virtualities (Qa, Qb smaller than a few GeV), the photon- 
photon cross section becomes dominated by soft interactions. Here one cannot use a per¬ 
turbative analysis. On the other hand, one may explore experimentally at what scales the 
breakdown of the perturbative result occurs, and how this is connected to the onset of the 
phenomenological “soft-pomeron” behavior. 

The theory that is available at present is leading logarithmic and therefore is affected 
by rather large uncertainties. These uncertainties can be parametrized in terms of two 
mass scales, the transverse scale that controls the running coupling and the longitudinal 
scale associated with the high energy logarithms. A next-to-leading logarithmic calculation 
would help determine these scales. Such a calculation could make the theoretical predictions 
much more precise. At the largest values of s, new effects related to unitarity and diffusion 
may become important. If so, an improved theory that deals with these effects would be 
testable at a future e + e _ collider. 
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APPENDIX A: THE BORN ORDER CALCULATION 

We start with the expression (|2.6|) for the amplitude corresponding to the graph in Fig. |2]. 
The overall charge factor in Eq. (|2.6| ) is 

el e 2 b e 4 Tr (t r t s ) Tr (t r t s ) = 32 a 2 a 2 s e 2 ) (2 tt) 4 , (Al) 

a,b \ q J 

where we have used the color trace Tr (t r t s ) = (l/2)5 rs . 

We use the mass shell constraints on the final quarks, pa and ps, to eliminate the integrals 
over their ” and “+” components, respectively, thus obtaining 
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d Z A j 2 . 


(2 7T 




-A = 


Pi 


2 *a gife 


(A2) 


d 4 p B 

(2vr) 4 


2vr^+(p|) 


2 


2 


p B (2 7T 


,-3 


% — 


P B 


2 sb gigs 


(A3) 


We use the mass shell constraints on the antiquarks to eliminate the integrals over the ’ 
and “+” components of the exchanged momentum k , as follows: 


/ YtF 2n6 +( ( y ( lA- Pa ~ k ) 2 ) 


2 (1 - z A ) qt 


+ ’ 


k ~ 


2gJ 


Qa + “ + 

i Z A 


P a , (Pa + k ) 5 


1-2,4 


(A4) 


/ YtT 2 ' tt6 +{( , 1b-Pb + * 0 2 ) 


2 Qb \ z b 1 ~ z b , 


2 (1 - 2b) g B 


(AS) 


Note that in Eqs. (|A|), (|A5|) we have neglected terms of order k + j\fs and k~/^/s with 
respect to unity, consistently with the high energy approximation. 

Eq. (|2.6|) can then be rewritten as 


\M\ 2 = 32 a 2 a 2 s ( 


(2 Y 


d 2 k d 2 p,4 d 2 pb 


x 


x 


■i 1 


d za 


,i i 


(2vr) 2 (2 7r) 2 (2 7r) 2 
d zb 11 


x 


Jo 2 7r 2 za (1 — za) Jo 2 7r 2 zb (1 — 2 b) 2 s ( k 2 ) 2 

Tr [ 7c (?*a + ft) 4a Ha zIaL ill ^ (4 a ~ ^a) 4] 

((pa + k) 2 + ie) ((p A - Qa) 2 -is) 

K 7“ (^b - #) 4b (4b - 4 b + f) i 3 (4 b - 4b) i 


Tr 


B 


«J>b - k) 2 + is) ((p B - 1 b) 2 - is) 


(A6) 


We now re-express the denominators and numerators of the amplitude (|2.6|) in the high 
energy limit. With the neglect of terms of order k + /-\/s and k~/^/~s, the denominators take 
the form 


(Pa + k) 2 ~ 


1 - 2,4 


(2,4 (1 — za) Qa + (pa + k)‘ 


(AT) 


(Pa ~ Qa ) 2 ^- (z A (1 - z A ) Qa + Pa) , 

^A 


(AS) 


and analogously 
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(Pb ~ k ) 2 ~ 


-to 


(to (1 — to) Qb + (Pb — k) 2 ) , 


(A9) 


(Pb - to ) 2 --(to (1 - to) Q% + Pi) 

to 


(A10) 


In the numerator, light-cone gluon polarizations are dominant at high energy, and therefore 
we are led to calculate the product of traces 


T+ + rr 


= Tr 
x Tr 


(All) 


i>A 7 + (4 + f) 4 a Ha - i> a -)) 7 + Ha ~ 4 a) 4 a 

i> B T 0 i> B - ft) 4 (4b ~1>B + t) 7 ' U> B ~ 4b) is 


The result reads 


T++T- b 


(A12) 


= 8 (q\f ( 424(1 - z a )za ■ Pa 2.4 ■ (pa + k) - pa ■ (pa + k)) 
x 8 (q B y- (4* fl (l - z b )e b ■ p B e B ■ (p B - k) - p B ■ (p B - k)) 


Substituting Eqs. (|A7|) - (|A10|) and (|A12|) in Eq. (|A6|), we obtain the expression (|2.7|) for 
the amplitude corresponding to the graph in Fig. ^]. 

The total cross section is arrived at by adding the contributions from the other graphs 
according to the replacements described in the text below Eq. O. and dividing by 2 s. 
The result reads 


= 128 « 2 « 2 


V 2 V f d2k d2 PA d 2 Pb /'' ,, f 1 

E e J (2 4 J (2 ( 2 7f)2 ( 2 7r ) 2 Jo d * A J 0 dzB (k 2 ) 2 


(A13) 


x 


X 


I [4to(l - to)to ■ P.4 to • (pa + k) - p 4 • (p A + k)] 

l [z A (1 - to) Qa + (PA + k) 2 ] [to (1 - to) Qa + Pa) 
[4to(l ^ to)(to • Pa) 2 - p\] \ 

[to (1 - to) Qa + Pa f J 
I [4to(l - to)to • Pb to ■ (pb - k) - p B • (p B - k)] 

l [z B (1 - to) Q% + (Pb - k) 2 ] [to (1 - to) Q% + Pb\ 
[4to(l - to)(to • Pb) 2 - p|] 1 
[z B (1 - to) Qb + Pb} 2 J 


This coincides with Eq. ( [2,1 1 ) in the text once the explicit expression ( 2 . 12 ) for G is used. 
Using Eqs. (|A 6 |) -( |ATT|) , the tensor Q ^ (k; Q 2 ) introduced in Eq. ( |3.6|) takes the form 


k; Q 2 ) = 


Ol Olq 


£• 


d 2 p 


4 (to ) 2 \q q JJ Jo 


d z 


(A14) 


' Tr [^ 7 + {j) + jt) Y {i - i> - ft) 7 + (j) - j) 7 "] 
x < —-rr^-—_ ___ r „ -—_ ___ + symm. 


[(p + k ) 2 + z (1 - z) Q 2 } [p 2 + z{l — z) Q 2 } 
where the additive symmetric terms are obtained from the replacements given in Sec. [IT]. 
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APPENDIX B: COMPARISON WITH THE CROSS SECTION FROM QUARK 

EXCHANGE 


The gluon exchange diagrams discussed in the text provide the dominant contribution 
to the photon-photon cross section in the high energy limit. They give rise to constant (in 
Born order) or logarithmic (in higher orders) terms at large s in the cross section. This 
appendix is concerned with quark exchange contributions, which vanish in the large energy 
limit. We examine quark exchange in order to estimate the energy at which gluon exchange 
becomes dominant. 


10 

5 


0 


-5 


10 

10 ° 10 1 10 2 10 3 

l/r 

FIG. 16. The r-dependence of the functions f\ and f 2 that enter the leading order expression 
for (j( q \ Eq. (|B4|). We report f± and /2 for two different values of Q\/Q 2 b . 

The leading-order term of quark-exchange type comes from the purely electromagnetic 
process 



- qI/qI = i 


- Qa./Qb = 5 

A-. - — 

l 

j _i_Li_i i i i I_i_i_i_i_i i i i I_i_i_i_i_i i i i 


1*(Qa) + 7 *(<? b ) -> q(p) + q(p) 


(Bl) 


This contribution is suppressed by a power of s at high energies, a ^ ~ 1/s. To express the 
cross section for this process, we parametrize the incoming photon momenta as in Eq. 
and introduce the variables 

c Q 2 a c Ql 

2 q\qB ’ ^ 2 q\q~ B 

in terms of which the total energy s has the expression 

s = 2 q a q B (i-£0(1-6?) • (B3) 
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FIG. 17. The T-dependence of the quark-exchange contribution to the 7*7* cross section in 
leading order. We take Qa = Qb, and we plot the rescaled cross section Qa Qb cd®)/ (« 2 (E q e$)). 
The dashed curve is the exact expression, Eq. f [iT~i|) , while the dot-dashed curve is the expression 
approximated for high energies, Eq. (B5). 


In the high energy region one has £4, -C 1, and s ~ 2 q^ as in Eq. 

The cross section for the process ( |B1| ) to order a 2 , averaged over the transverse photon 
polarizations, has the form 


r(9) _ 27Ta 2 


cr™' = 


VUJb 


QaQb 2(1 -£a£b) 


[/1 + /2 In (1/(6^))] 


(B4) 


where f\ and f-z are rational functions of £a, £b, and are plotted in Fig. 0 versus the variable 
r = \/£aXr for different values of the ratio p = £a/£b- 

The r-dependence of the cross section ( |B4| ) is reported in Fig. 13 for the case of equal 


virtualities. The cross section vanishes at the kinematic threshold r = 1, it has a maximum 
around r ~ 10 _1 , then it falls off and vanishes for r —> 0 (corresponding to high energy) like 
r In r, according to the asymptotic formula 


a ( Q ) ^ 27ra 4r (ln(1/r) _ 1} 


Qa Q 


rCl 


(B5) 


B 


The power suppression with r at small r is the one expected from the exchange of a spin- 
1/2 line in the t-channcl. The logarithmic enhancement is associated with the integration 
over the region of small angles at the splitting vertex 7* —> q q in the limit of small photon 
virtuality. The behavior of the cross section is qualitatively the same in the case of unequal 
virtualities. 


42 








REFERENCES 


[ 1 ] 


[ 2 ] 


[3] 

[4] 

[5] 


[ 6 ] 

[7] 


[ 8 ] 

[9] 

[ 10 ] 

[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 

[ 22 ] 

[ 23 ] 


L.N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E.A. Kuraev, L.N. Lipatov and V.S. 
Fadin, Sov. Phys. JETP 45, 199 (1977) ; I. Balitskii and L.N. Lipatov, Sov. J. Nucl. 
Phys. 28, 822 (1978). 

H. Abramowicz, plenary talk at ICHEP96 (Warsaw, Jnly 1996), in Proceedings of the 
XXVIII International Conference on High Energy Physics, eds. Z. Ajduk and A.K. 
Wroblewski, World Scientific, p.53. 

DO Collaboration, Phys. Rev. Lett. 77, 595 (1996). 

S.J. Brodsky, talk at Workshop on High Energy e + e _ Colliders, Brookhaven National 
Laboratory, May 1996. 

F. Hautmann, talk at ICHEP96 (Warsaw, July 1996), preprint OITS 613/96, in Pro¬ 
ceedings of the XXVIII International Conference on High Energy Physics, eds. Z. Ajduk 
and A.K. Wroblewski, World Scientific, p.705. 

S.J. Brodsky, F. Hautmann and D.E. Soper, Phys. Rev. Lett. 78, 803 (1997). 

P. Aurenche, G.A. Schuler et ah, Report on “77 Physics” in Proceedings of the Work¬ 
shop “Physics at LEP2”, eds. G. Altarelli, T. Sjostrand and F. Zwirner, CERN 1996-01, 
Vol.l, p.291. 

A.H. Mueller, Nucl. Phys. B415, 373 (1994); A.H. Mueller and B. Patel, ibid. B425, 
471 (1994). 

I. Balitskii, Nucl. Phys. B463, 99 (1996). 

J. Bartels, A. De Roeck and H. Lotter, Phys. Lett. B 389, 742 (1996). 

V.M. Budnev, I.F. Ginzburg, G.V. Meledin and V.G. Serbo, Phys. Rep. 15 C, 181 
(1975). 

F. E. Low, Phys. Rev. D 12, 163 (1975); S. Nussinov, Phys. Rev. Lett. 34, 1286 (1975), 
Phys. Rev. D 14, 246 (1976); J.F. Gunion and D.E. Soper, ibid. 15, 2617(1977). 

J.D. Bjorken, J. Kogut and D.E. Soper, Phys. Rev. D 3, 1382 (1971). 

J.D. Bjorken and J. Kogut, Phys. Rev. D 8, 1341 (1973). 

L.L. Frankfurt and M. Strikman, Phys. Rep. 160, 235 (1988). 

S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B 242, 97 (1990), Nucl. Phys. 
B366, 135 (1991); J.C. Collins and R.K. Ellis, ibid. B360, 3 (1991). 

G. P. Lepage and P.B. Mackenzie, Phys. Rev. D 48, 2250 (1993); S.J. Brodsky, G.P. 
Lepage and P.B. Mackenzie, ibid. 28, 228 (1983). 

J.D. Bjorken, preprint SLAC-PUB-7341, presented at Snowmass 1996 Summer Study 
on New Directions for High Energy Physics, e-print archive [hcp-ph/9610516 . 

J. Bartels and H. Lotter, Phys. Lett. B 309, 400 (1993). 

A.H. Mueller, Nucl. Phys. B437, 107 (1995). 

G.A. Schuler and T. Sjostrand, Zeit. Phys. C 68, 607 (1995), Phys. Lett. B 376, 193 
(1996); M. Gliick, E. Reya and M. Stratmann, Phys. Rev. D 51, 3220 (1995). 

P.D.B. Collins, An introduction to Regge theory and high energy physics , Cambridge 
University Press, Cambridge, 1977. 

OPAL Collaboration, contributed paper pa03-007 at ICHEP96 (Warsaw, July 1996); 
J.A. Lauber (OPAL), in Proceedings of the XXVIII International Conference on High 
Energy Physics, eds. Z. Ajduk and A.K. Wroblewski, World Scientific, p.725. 


43 



[24] S. Kuhlman et al., Physics and Technology of the Next Linear Collider, Snowmass 1996 
Report, e-print archive [hep-ex/9605011 . 


44 





